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We investigate the direct detection phenomenology of a class of dark matter (DM) models in 
which DM does not directly interact with nuclei, but rather the products of its annihilation do. 
When these annihilation products are very light compared to the DM mass, the scattering in direct 
detection experiments is controlled by relativistic kinematics. This results in a distinctive recoil 
spectrum, a non-standard and or even absent annual modulation, and the ability to probe DM 
masses as low as a ~ 10 MeV. We use current LUX data to show that experimental sensitivity to 
thermal relic annihilation cross sections has already been reached in a class of models. Moreover, the 
compatibility of dark matter direct detection experiments can be compared directly in i?min space 
without making assumptions about DM astrophysics, mass, or scattering form factors. Lastly, when 
DM has direct couplings to nuclei, the limit from annihilation to relativistic particles in the Sun can 
be stronger than that of conventional non-relativistic direct detection by more than three orders of 
magnitude for masses in a 2-7 GeV window. 


Introduction - While very little is known about 
Dark Matter (DM), its cosmological abundance is exper¬ 
imentally quite well-determined: UtcoMh"^ = 0.1199 ± 
0.0027 [T]. An appealing framework for understanding 
the relic abundance of Dark Matter (DM) is thermal 
freeze-out [2]- Number-changing interactions in the early 
universe, XX o (SM)SM keep DM in thermal equilib¬ 
rium with the SM bath, until the rate of these annihila¬ 
tion processes drops below the rate of Hubble expansion. 
After this point the abundance of DM is essentially fixed 
at, IlcDMh'^ ^ 0.12(6 X 10“^® Cm^ S“V(o'annWez)), 
singling out a characteristic annihilation cross section 
(o’annVrei) for thermally produced DM to yield the ob¬ 
served abundance. This scenario is attractive in that 
it provides a simple and elegant framework for the relic 
abundance that can be tested in a variety of ways, in¬ 
cluding direct detection (DD) [3]. However, current con¬ 
straints from DD rule out many of the simplest models 
of thermal relic DM, which may indicate a modification 
of the above picture. 

In this paper we investigate a modification of thermal 
DM which alleviates the tension between DD constraints 
and the thermal relic hypothesis, while making unique 
predictions for DD. In particular, we take the abundance 
of DM, X, to be determined by the annihilation pro¬ 
cess A A -O' YY, where A is a much lighter dark sector 
species. The interactions of the dark sector state Y with 
ordinary nuclei allows for a unique test of the scenario at 
DD experiments. The resulting DD phenomenology of 
this class of models is distinctive, owing to the fact that 
(1) the scattering partner of the nucleus is relativistic, 
rendering the kinematics of scattering completely differ¬ 
ent and (2) it is the flux of the scattering partner Y that 
determines the rate of events at a detector rather than A. 
Both of these features have novel consequences not con¬ 
sidered in the literature of “model- independent” direct 
detection analyses nig. 

As loop processes will always engender scattering of A 
on nuclei at DD, we will focus on DM masses less than 


~GeV such that the non-relativistic scattering of A does 
not produce detectable nuclear recoils above a detector’s 
O(keV) threshold. Similar scenarios have recently been 
investigated in [SH with a focus on the Cherenkov sig¬ 
nals at Super-Kamiokande and IceCube. 

In this paper we employ current LUX nni limits to 
demonstrate that DD experiments are sensitive to ther¬ 
mal relic annihilation cross sections for galactic center an¬ 
nihilation of DM in a window of DM masses from 10 MeV 
to 1 GeV. Direct detection has historically been muddled 
by multiple conflicting data sets. To combat this, we il¬ 
lustrate how current and future direct detection data can 
be easily analyzed for compatibility in this framework by 
mapping results to Amin-space. Additionally, we investi¬ 
gate the testability of such relativistic scattering models 
where the signal is dominated by DM accretion and an¬ 
nihilation within the sun. 

Annihilating DM in the Galactic Center - Two 
potential sources of DM annihilation are annihilation 
from the Galactic Center and annihilation within the 
Sun. A key difference between these two is that the lat¬ 
ter relies on a stable balance between the accretion and 
evaporation rates of DM interacting with nucleons inside 
the Sun. We first consider galactic center annihilation 
since this does not require a build up of DM in the Sun 
and hence requires fewer assumptions. 

For simplicity consider “2-to-2” annihilation, AA —> 
YY. Then the differential rate (per unit detector mass) 
at a direct detection experiment is, 

= — r dEy^(^), (1) 

dEfi rriN dEy V dE^ J 

where ‘by is the local flux of A’s, Ej^iniER) = 
•y/ mNEji/2 is the minimum energy to produce a recoil 
of energy Er, and ^ = 25 {Ey -mx)- 

For simplicity, we adopt a contact interaction between 
Y and a quark q of the form Oqy = Gy {Y^^Y) {qj^q) 
where Gy is the effective coupling. By analogy with 
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FIG. 1. : LUX Limits on XX —>■ Fsmi^sm, for which 
Gy ~ Gf = 1.2 X 10~® GeV“^, and a model in which DM 
annihilates to relativistic pairs YY which scatter on nuclei, 
YN YN, via Eq. In the shaded gray band we vary Gy 
in the interval (7 x 10^ — 5 x 10^) Gf- Propagation effects 
have been accounted for in determining the flux of Y at the 
underground site of the LUX detector (see Appendix). 


neutrino-nucleus elastic [^scattering the differential cross 
section is m, 


dcjYN 

dEn 


^A'^mNF^iER) 



( 2 ) 


where F{Eji) is Helm’s nuclear form factor [T^], and the 
coefficient is for the coherent enhancement of scatter¬ 
ing with equal rates on protons and neutrons. 

Next we must determine the local flux of Y. The 
flux of Y particles from DM annihilation in the 
Galactic Center is estimated as [7], = 1-6 x 

in— 2„ —2„ — 1 f (^xx-»yy Vr-ed ^ { 20 MeV A 

iU cm S (^5x10-26 cm3s-i ^ nix ) ' 

Given this flux, and a model of T-nucleus interactions, 
the only remaining parameter to determine is the an¬ 
nihilation cross section, which we take to be a free pa¬ 
rameter, to be determined from data. First let us take 
a minimal choice by relying on the SM to furnish the 
interactions of Y with the nucleus. This immediately 
singles out the neutrinos as the only SM possibility for 


^ Note that in m it was estimated that the inelastic scattering 
cross section is small at low-energies compared to the elastic 
cross section, cr(elastic)/c7{inelastic) ~ /(Ey, where the 

nuclear radius is, Rn ~ (10 MeV)“^. Thus for a Xenon tar¬ 
get nucleus, inelastic scattering is sub-dominant for Ey < GeV. 
Given our focus on sub-GeV DM we will ignore inelastic pro¬ 
cesses in this paper. 


Y. The elastic, spin-independent scattering of SM neu¬ 
trinos with nuclei can be computed using Eq. with the 
replacement, Gf{N/2)‘^, where N is the num¬ 

ber of neutrons and Gf is the Fermi constant. We see 
in Fig. that with present LUX data, the resulting sen¬ 
sitivity to the annihilation cross section is weak, being 
orders of magnitude away from thermal relic sensitivity. 

On the other hand, DM could well annihilate to non- 
SM particles that have larger than electroweak-size inter¬ 
actions. Two generic classes of models serve as examples: 
models of gauged baryon number [T3H13 and so-called 
“Higgs portal” models . Gauged baryon number is 

motivated by the stability of the proton, which in the SM 
remains a mystery and may indicate that baryon num¬ 
ber is in fact a gauge symmetry. This is one of the few 
phenomenologically viable “portals” connecting the dark 
and visible sectors, as it does violate any of the approxi¬ 
mate symmetries of the SM. In addition, the Higgs por¬ 
tal, .if D |(/)p|iFp (where H is the SM Higgs and </> a dark 
sector scalar), represents a rather generic possibility for 
connecting the dark and visible sectors. In this case it is 
natural for scattering on nuclei to be enhanced relative 
to electrons since the couplings scale with the SM-Higgs 
Yukawa couplings. 

For illustration, we can make use of a simplified model 
for quark-Y interactions via the exchange of a light vector 
or scalar </>: 

.ifv D [gqqitiq + gYYj^Y) , (3) 

.ifs D ^ {gqqq -b , (4) 

where gY-,gq are the couplings of SM quarks and Y to 
the mediator. In terms of these couplings, the effective 
parameter is, Gy = (gqgY) and the relevant con¬ 
straints for the two models are discussed respectively in 
the Appendix. 

As a benchmark we take Gy in Fig. [l] to vary in the 
interval (7 x 10^ — 5 x 10^) Gp- We highlight that the 
values of Gy are well within the constraints allowed by 
“missing energy” collider limits [33 [H] ■ Larger Gy are 
permitted by collider limits, though the flux of Y parti¬ 
cles becomes strongly suppressed (see the Appendix for 
a discussion of this effect). Nonetheless, we see in Fig. 
that models of this type are already being probed by di¬ 
rect detection and can in particular exclude thermal relics 
in the 10 MeV - 0.5 GeV window, yielding a novel probe 
of thermal DM. Future constraints will cut further into 
thermal relic territory. 

It is important to observe that annihilation of DM 
to relativistic states from the galactic center predicts 
no sizeable annual modulation. In the case of annihi¬ 
lation from the Sun however, the annual modulation is 
known to peak in January due to the eccentricity of the 
Earth’s orbit. Thus solar neutrino signals in direct de¬ 
tection experiments predict a nearly maximally “wrong” 
phasej^ with respect to the expectation from light non- 
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We note that the phase can be reversed for heavy DM at a high 
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relativistic DM of June 2""^ [3^. This expectation can be 
violated however when the annihilation product Y expe¬ 
riences flavor oscillations on C>(AU) length scales as in 
for example HD 131133 though this requires very small 
mass-splittings, Am^ ^ 10“^° eV^. 

It is important to observe that the ability to scatter on 
nuclei does not induce any physics which would allow Y 
to decay. Given the generic stability of both X and Y, we 
must be sure that their total abundance does not exceed 
the observed value, flcDAih^ — 0.2. For simplicity, we 
will work in the limit that Y forms a sub-dominant com¬ 
ponent, i.e. fly ^ fix — flcDM- This can be naturally 
arranged when Y is similar to a neutrino and freezes-out 
when it is relativistic, i.e. my < eV. Non-relativistic 
freeze-out of Y can also lead to a small relic abundance 
when its annihilation cross section is large [ 3 . 

Another potential constraint on these models is the 
additional radiation energy density they generate during 
the Big Bang, parameterized by, ANes = ^ = 
where the photon and neutrino temperatures are re¬ 
lated by, Ty = (4/11)^^^ Ty. Due to the annihilation 
of degrees of freedom from the standard model plasma, 
the temperature of the dark sector relative to the stan¬ 
dard model sector is suppressed via dilution, Tdark = 
( 5 *. post-BBN/ 5 *, dark freeze where the number 

of standard model degrees of freedom (DOFs) after Big 
Bang Nucleosynthesis (BBN) is post-BBN = 3.36, the 
number of standard model DOFs after dark freeze out is 
9*, dark freeze outj and Ty is the photon temperature. As¬ 
suming Y is the lightest stable particle in the dark sec¬ 
tor, so that heavier dark sector DOFs re-heat the bath 
of Y radiation as they annihilate. Under the assumption 
that entropy is conserved, gosT^ark — dvTy, where g^s 
count the total DOFs in the dark sector. Typical values 
of AA^eff — 0.2 for dark freeze out Corresponding to the 
QCD phase transition and minimal additional degrees of 
freedom in the dark sector (only X and Y). This is well 
within the allowed constraints on AAeff [DI3H] • 

Direct Detection in space - Direct detection 

involves a unique combination of particle physics, nuclear 
physics, and astrophysics. The kinematics of scattering 
in the non-relativistic case are controlled by the minimum 
DM particle velocity, Vrain{ER), required to produce a nu¬ 
clear recoil of energy Er . In the absence of unknown form 
factors, all experimental data can be mapped into Umin- 
space at each DM mass and compared without specifying 
the nature of the astrophysical distribution or density of 
DM [331 HO]- These “halo-independent” methods have re¬ 
ceived signiflcant attention |1TH55] . We generalize these 
methods to cover relativistic scattering as well, where the 
“halo-independence” here comes from the absence of spe¬ 
cific assumptions regarding the local DM density, density 
profile, velocity distribution, and annihilation source. 

As can be seen from Eq. Q and the form of Amm(£'/?), 
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FIG. 2. DM mass independent comparison of direct detec¬ 
tion data under the assumption of relativistic scattering. Here 
we include DAMA modulation amplitude from [53 and the 3 
CDMS-Si events [58] . along with constraints from LUX m 
and SuperCDMS |59| . 


the relativistic scattering case allows a comparison of 
data which is independent of the DM mass. In the case 
of a claimed detection, using Eq. Q we can divide out 
the nucleus-specific quantities 

g{E^in) = 2gl{A^F\ER))-^dR/dER, (5) 

to immediately obtain the preferred g range in Fig. [^ 
Finally, since the integrand in Eq. [^ is strictly positive 
we can derive conservative limits on g as in m by as¬ 
suming a step function form for g{E^i^). One can view 
this procedure as mapping direct detection rates to the 
{g — Ejnin) plane, which we refer to as “Emin-space” for 
brevity. The form of Eynin{ER) has the interesting effect 
of strongly suppressing the sensitivity of experiments em¬ 
ploying heavy target nuclei. It is also interesting to ob¬ 
serve that LUX uni and a relativistic DM interpretation 
of DAMA [33 and CDMS-Si [33[ data are fully compat¬ 
ible, though essentially ruled out by the recent Super¬ 
CDMS data [33]. Clearly, allowing for isospin-violation 
in order to suppress the sensitivity from Germanium¬ 
scattering would result in the positive signals seen by 
CDMS-Si and DAMA and the null results of LUX and 
SuperCDMS to be compatible. 

Let us pause to highlight the relevance and general¬ 
ity of the halo-independent method employed here. In 
contrast with non-relativistic scattering, here the veloc¬ 
ity distribution matters very little for the rate of events. 
However, now the astrophysical uncertainty is more fun¬ 
damental in the sense that the souree of the flux is un¬ 
known, i.e. the Galactic Center, the solar interior, etc. 
Moreover, even after specifying a source there exist large 
uncertainties in the spatial distribution. This method is 
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independent of these sizeable uncertainties. Finally, in 
addition to being DM mass independent, this method is 
also automatically independent of the energy dependence 
of the F-nucleus scattering and the spectrum of the Y 
particles. The generality of this method is illustrated by 
the model of “baryonic neutrinos” which was proposed to 
account for DAMA’s annual modulation mi 131133 iS]- 
In this case, despite the fact that the source of the rel¬ 
ativistic Y’s is completely disconnected from DM, the 
Aniin-space representation is valid and allows for a com¬ 
plete comparison of experiments as in Fig.[^ 

DM Annihilation from the Sun - For solar annihi¬ 
lation to dominate over the contribution from the Galac¬ 
tic Center, the Sun must contain a large quantity of cap¬ 
tured DM. In a symmetric DM context, the solar DM 
abundance has the time evolution Nx = Cxn — 
such that solar annihilation fluxes are roughly <I >0 ~ 

, where Ca is the annihilation rate, Nx is the num- 
ber of DM particles in the Sun, and Rau is the Earth-Sun 
distance. Assuming that only annihilation and nuclear 
capture play a role we can specify a model of scatter¬ 
ing of the Y states on nuclei and then derive bounds on 
the DM-nucleus cross section. This assumption is valid 
in the regime where evaporation of DM out of the sun is 
negligible, N^q — \/CxnICa-, such that the flux depends 
only on the capture rate. 

We again adopt the model of DM annihilating to dark 
sector Y states that interact with quarks via Eq. Q and 
take gv = 0.1 and my = 50 MeV. This yields the result 
shown in Fig. At low DM mass the limit cannot be 
trusted, as sufficiently light DM is prone to evaporation 
from collisions with solar nuclei [5TH53] . We note that in 
models where DM experiences significant self-interaction 



FIG. 3. Probing DM-proton interactions from the annihila¬ 
tion in the Sun. For reference we show the cross section limits 
from LUX and SuperCDMS derived under the assumption of 
ordinary non-realistic scattering. Here the scattering of the 
relativistic state is via vector with Gy = 7 x 10^ Gf. 


the abundance of DM in the Sun can be much larger, 
which can strengthen the limit in Fig. [^significantly |64F 
inzi- We also leave for future work the extension of the 
framework considered here to an asymmetric DM sce¬ 
nario (see e.g. [SSHTO] and [55)1. 

Discussion and Summary- In summary, this 
work has investigated the sensitivity of direct detec¬ 
tion searches to dark matter annihilation. Thermal relic 
dark matter sets a natural scale for the thermally aver¬ 
aged DM annihilation cross section around {(JannVrei) — 
6 X 10“^® cm^ s“^. This scale can be searched for in 
CMB [71], gamma-ray |72|, and even neutrino data [73] . 
Both CMB and gamma-ray data have breached ther¬ 
mal relic sensitivity for light DM masses. Though these 
constraints have sizable astrophysical uncertainties, they 
may indicate that light DM requires non-SM modes of 
annihilation. Here we have studied models in which DM 
annihilates to a light, non-SM state that can scatter elas¬ 
tically on nuclei and deposit a detectable recoil energy. 
Models of the type considered here retain the appeal of 
the thermal relic hypothesis while remaining experimen¬ 
tally verifiable. We have furthermore demonstrated that 
in this class of indirect annihilation searches, all astro- 
physical uncertainties can be “integrated out ” m and 
experimental sensitivities can be directly compared. 

This work could be extended to include electronic scat¬ 
tering at direct detection m, though the reduction in the 
Cherenkov threshold for electrons implies that Super-K 
limits extend to much lower masses for leptophilic mod¬ 
els [7]. The most similar studies to our own which have 
been recently carried out assumed that DM interacts 
with the SM through a kinetically mixed photon, im¬ 
plying both hadronic and electronic couplings [UIH] . In 
this case, large volume detectors like Super-Kamiokande 
and IceCube yield very strong limits. In contrast, we 
are interested in a complementary portion of the param¬ 
eter space compared to Hi in that we have focused on 
hadronic models where: ( 1 ) the annihilation products are 
nearly massless compared to nuclear recoil energies and 
(2) light DM masses which are near or below Cherenkov 
threshold and thus difficult to probe at Super-K. 

Lastly, we stress that the “effective field theory” of 
DD proposed in [ 1 [S] does not encapsulate the scenario 
described in this paper and should be extended to include 
generalized relativistic scattering. 
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APPENDIX 


A. Stopping Effects 


The mean free path of the Y particles as they enter the 
Earth is, A = (na)~^. The total T-proton cross section 
is 


ay-p ~ 4 X 10 cm^ 


/ Sr Yf Gy 

\10Mev) ^vTxIO^Gf 


2 


( 6 ) 

For a conservative estimate, we take the Earth to be en¬ 
tirely composed of Iron and obtain njv = p^/mpi, ~ 
5.5 X 10^^ cm“^ where = 5.5 g cm“^. Thus the mean 
free path of Y is 


A = -~ 7000 km 

riNayN 


/10 MeV 

V Ey 


2 


^7 X lO^Gi.^ 


(7) 

This can be approximately incorporated into Fig.l by 
replacing the flux with —>• $vx (exp (—£(0)/A)) where 
the brackets indicate averaging over arrival directions and 
£{9) is the chord length an incoming Y particle traverses 
in the Earth prior to arrival at the detector as a function 
of the incident angle, 9. 


B. Higgs portal Example 


In a Higgs portal implementation, we have .jSf D 
b\(j)\'^\H\'^ + g(j)YY. Under the assumption that the scalar 
potential breaks the global 1/(1), it will develop a vacuum 
expectation value and mass mix with the SM Higgs, H. 
The precise value of the mixing angle depends on the 
scalar sector, which is a model-dependent feature. As is 
well-known, the Higgs portal is constrained by the invis¬ 
ible branching ratio of the Higgs, which is constrained to 
be < 26% [31] (though see |7S] for weaker, direct con¬ 
straints). In our case, the Higgs has two new contribu¬ 
tions to its invisible width, h —> (jxp and h —> TT, with 
widths 


r{h —>■ (j)(f>) 

T{h YY) 



( 8 ) 

(9) 


We find that these two processes respectively require, b < 
7.1 X 10-3 and gsm9 < 1.7 x IQ-^. 

For sufficiently large (j) mass, we can parameterize 
the Y-N interaction as contact and write the vertex as 
Gy(NN)(YY) where 


Gy 


[g sin9)fN 


7x 103 Gf 


fgsin9\ /0.2 GeVV 

I J 


( 10 ) 

where fjv = 0.345 [32] is the effective coupling of the (j) 
particles to nucleons. We conclude that the Higgs portal 
is a viable possibility for the interactions of Y and nuclei. 
Note further that the coupling to leptons is naturally 
Yukawa suppressed. 


C. Gauged Baryon Number Example 


Interactions of the type described by arise in 
models of gauged baryon number [niiMi 137], gauged 
U{1)b-l [Sni [ZS], and kinetic mixing |60| . though it is 
only in the first case that couplings to leptons are ab¬ 
sent. A detailed exploration of the allowed parameter 
space of Eq. Q is beyond the scope of this work, but 
see m for a useful compendium of constraints. 

In the U{1)b case the dominant coupling to hadrons is 
explicitly guaranteed by the baryon symmetry. The La- 
grangian includes J^’b D + grYj^iY). Light 

baryonic gauge bosons are constrained by the invisible 
decay width of the T(1S) [TH] US], which is constrained 
to have a branching ratio, H7?. (T(IS)—>■ invisible) < 
3x 10“^ [77], missing energy from the LHC/Tevatron [33] 
and rare radiative decays of the light mesons [22] . As an 
example we take my < 100 MeV, where we hnd the 
mono-jet constraints from CDF to be strongest. In the 
limit my my, these searches constrain the combina¬ 
tion of parameters gq^BRiVp —)■ YY) < 0.02. 

Mapping again to the effective parameter Gy appear¬ 
ing in Eq.([^ we find 


Gy 


» cs 7 X 103 
m^ 



/35 MeV 
V my 


, ( 11 ) 


and thereby conclude the “baryonic portal” is another 
viable model yielding the phenomenology outlined in the 
main text of the paper. 
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